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Baryons in the large TV limit of (f +f )-dimensional Gross-Neveu models with either discrete or 
continuous chiral symmetry have long been known. We generalize their construction to the case 
where the symmetry is explicitly broken by a bare mass term in the Lagrangian. In the discrete 
symmetry case, the exact solution is found for arbitrary bare fermion mass, using the Hartree-Fock 
approach. It is mathematically closely related to polarons and bipolarons in conducting polymers. 
In the continuous symmetry case, a derivative expansion allows us to rederive a formerly proposed 
Skyrme-type model and to compute systematically corrections to the leading order description based 
on an effective sine-Gordon theory. 
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I. INTRODUCTION 

The original Gross-Neveu models Q| are asymptoti- 
cally free, self-interacting fermionic field theories in 
dimensions with Lagrangian 

£ = ^ (i7'^5^ - mo) V- + ^5' + ??(Vii7V)'] • (1) 

We suppress the flavor indices {i = \...N) and shall only 
consider the 't Hooft limit N — > 00, Ng^ ~ const. The 
model with discrete chiral symmetry [ry = in Eq. iQ] 
will be referred to as Gross-Neveu model (GN), the one 
with continuous chiral symmetry (ry = 1) as Nambu- 
Jona-Lasinio model in two dimensions (NJL2) through- 
out this paper. Apart from being of theoretical interest 
in their own right, these models have been useful as a the- 
oretical laboratory for new methods and algorithms and 
to describe quasi-one-dimensional condensed matter sys- 
tems. Whereas the massless models (toq = 0) have been 
rather comprehensively studied by now (see l^jQ for per- 
tinent review articles and ^ ^ ^ for a recent update on 
the GN model), results about the massive versions with 
an explicit symmetry breaking term (toq 7^ 0) are scarce 
and somewhat scattered through the literature. In view 
of the fact that in nature quarks are massive and a lot 
of effort is presently devoted to computing chiral correc- 
tions, we believe that it is both worthwhile and timely 
to study more systematically the massive Gross-Neveu 
models. In order to motivate our specific work, let us 
briefly summarize the present status of this field. 

The gap equation which determines the vacuum and 
the dynamical fermion mass m can easily be generalized 
to finite bare quark mass mg. Apart from to, a sec- 
ond physical, renormalization group invariant parameter 
appears, TOo/iVg^. It enters the renormalized vacuum 
energy and observables like the qq scattering amplitude 
or the mass of the "pion" in the NJL2 model The 
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CT-meson is unbound and disappears from the spectrum 
for any toq 7^ 0, the 7r-meson becomes massive in the fa- 
miliar way (Gell-Mann, Oakes, Renner relation). If one 
assumes unbroken translational invariance, one can also 
infer the phase diagram in the (r,/x) plane % However, 
since we now know that this assumption is not justified 
at Too — 6|, this analysis cannot be trusted and needs 
to be repeated. 

As far as baryons in the massive Gross-Neveu mod- 
els are concerned, we are only aware of two variational 
calculations so far. Salcedo et al. Q have studied both 
the 't Hooft model (large N QCD2 [13) and the NJL2 
model, comparing a lattice Hartree-Fock calculation with 
a variational ansatz in which the chiral phase of the 
fermions serves as effective low energy field. In the limit 
of small but finite bare quark masses, the authors find 
a sine-Gordon theory, the two-dimensional analogue of 
the Skyrme model where baryon number is generated 
through a topologically non-trivial pion field configura- 
tion. This elegant approach yields the correct depen- 
dence of the pion mass on the quark mass in a very simple 
manner and can readily be generalized to finite temper- 
ature and chemical potential 2, ffj. In the case of the 
't Hooft model, the authors of Ref. (9| compare the soli- 
ton approach with the full Hartree-Fock calculation and 
notice that it works very well, even away from the chiral 
limit. Unfortunately, such a comparison was not done for 
the NJL2 model, so that the accuracy of the sine-Gordon 
approach cannot be assessed in the case at hand. 

In the non-chiral GN model, Feinberg and Zee have 
performed a variational calculation, based on the scalar 
potential of the standard (toq = 0) kink-antikink baryon 
[13] . They compute the energy of the baryon and discuss 
the limiting cases of small and large bare quark masses. 
They conclude that their ansatz does not satisfy the sad- 
dle point equation. Claims made in this paper about the 
non-existence of static bags in the massive GN model 
have later been retracted by the authors [l^ . 

Since one does not know the accuracy of these varia- 
tional calculations, they are insufficient if one wishes to 
use the Gross-Neveu models as testing ground for new 
theoretical approaches. Moreover, as should be clear 
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from the above short summary, the general situation con- 
cerning the massive Gross-Neveu models is far from sat- 
isfactory. We therefore decided to reconsider this whole 
issue more systematically. In the present work, we begin 
such a study with a fresh look at the baryon in both the 
massive GN and the NJL2 models. Obviously, this is a 
prerequisite for a future study of the full phase structure 
of these models. In Sect. II, we will present the analyti- 
cal solution for the baryon in the discrete chiral massive 
Gross-Neveu model for arbitrary quark masses, obtained 
along lines very similar to Ref. |l2l |. In Sect. Ill, we turn 
to the NJL2 model and show that the Skyrme-type ap- 
proach 1^ can be identified with the leading order term of 
a systematic derivative expansion. We will also compute 
the three following higher order corrections in the deriva- 
tive expansion in closed analytical form. It seems to us 
that the problem of baryonic matter at finite tempera- 
ture away from the chiral limit should be tractable owing 
to these new insights, but we leave such an investigation 
for the future. Sect. IV contains a brief summary and 
our conclusions as well as a comment on the quantitative 
relationship between baryons in the massive GN model 
on the one hand and polarons, bipolarons and excitons 
in conducting polymers on the other hand. 



II. BARYONS IN THE MASSIVE GN MODEL 
— EXACT RESULTS 



Here, we basically repeat the variational calculation 
of Ref. 0, using the Hartree-Fock approach. We sim- 
ply write down a scalar trial potential and prove its self- 
consistency. Since our trial potential has the same shape 
as the mo = baryon with partial filling of the valence 
level (a kink-antikink potential well _13J), we can take 
over many results from the corresponding Hartree-Fock 
calculation in the chiral limit TE] . We refer the reader to 
this last paper for more details and will concentrate on 
the differences caused by the bare quark mass. 

The trial scalar potential reads 



S{x) — ni[l + y (tanh^_ — tanh^+)] 



(2) 



with 



ymx ± — arctanhy. 



where m is the physical fermion mass in the vacuum, 
y G [0, 1] the only variational parameter. The potential 
is the same as in the mo = case, except that we will 
get a different relation between the parameter y and the 
occupation fraction v = n/N of the valence state. We 
first evaluate the baryon mass relative to the vacuum. 
The calculation of the sum over single particles is identi- 
cal to that in Ref. 15]. Hence the contribution from the 



negative energy continuum states (Eik) = + m^ ), 
2Nym f^'"^ dk 1 



A/2 2^ E{k) 

2Nm I J\ - y2 

H V 1 ~ y arctan — 

TT y 

2Nym f^^-^ A 



+ 



arctan 



and the discrete states 



ASdiscr = -N{1 - iy)my/l-y^ 



(4) 



(5) 



can be taken over literally. Only the double counting 
correction of the interaction energy characteristic for the 
Hartree-Fock approach gets modified due to the bare 
quark mass mo. 



AE^d.c. = / dx 



{S - ?7io)^ - (m - mo)^ 



25^ 



= ^Jdx{{S^-m^)-2mo{S-m)}.{6) 

Moreover, when eliminating the bare coupling constant, 
we must now use the finite mo gap equation, 
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(7) 



where we have dropped mo-terms not matched by 
ln(A/m). This reproduces the previous result plus a fi- 
nite correction, 



^E^ ^ ^ I dx{S'^ - rr?) 

2tt m 



2iT ra m J 



(8) 



Carrying out the integrations and introducing the phys- 
ical parameter 



mo , A 

7 = m — 

7rm m 



(9) 



(3) then yields 



2NvTn A 

AEdc = ^In 2A^m7(?/- arctanhy). (10) 

TT m 

Collecting the results Q), lO and (fTn|l . we get the fol- 
lowing answer for the (variational) baryon mass 



^ . ^ + ^yr^^arctan^iHZ (n) 

N TT TT y 



-(1 — v)ra\J 1 — — 2mj{y — arctanh y). 
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We now choose v (i.e., the fermion number) and vary Mb 
with respect to y 



dy 



= 



(12) 



or (discarding the trivial solution y ^ 0) 



V 1 - 2/ - arctan — = -/y. (13) 

\^7r y 2 J 

Introducing the angle 9 via y — sin 9 (0 < 9 < 7r/2), we 
obtain 



— = — h 7 tan ( 

2 TT 



(14) 



If we eliminate i' from Eq. 111|) with the help of Eq. (|14|l . 

the baryon mass at the minimum finally becomes 

^^^^ = sin0 + 2m7arctanh(sin^?). (15) 

N TT 

The last two equations agree with Ref . . 

Now consider the self-consistency condition for the 
condensate and scalar potential in the form 



{S - Too) 



(16) 



The r.h.s. of Eq. (|16|l gets contributions from the discrete 
states 

occ 

Y.i'i' = (1 - i^)y Vl - (tanhC- - tanhe+) 



discr 



= (1-^)^^V^(5~to) 

and from the negative energy continuum 
^^/^ dk 1 



(17) 



m^y{l-y^) (18) 
1 



A/2 2^ E{k) 

(*^"^^--'""^^+)/ii?(fc)(P+TOV) 



^ -^lnA_(^VrZZ,rctan^i3. 

TT TO TT y y 

where we have again taken over results from Ref. |l5l |. 
The l.h.s. of Eq. (|16|l can be rewritten with the help of 
the gap equation 0, 

(S — toq) S" a S — m toq a 
= m m — 



TT TO TT TO TO 

— In ~ m)"f. 

TT TO 



(19) 



Combining Eqs. H17|l . (|18() and H19() . the self-consistency 
condition assumes the form 



■ arctan -^^ h 7 = 



2 y TT y 



y 



(20) 



which coincides with the variational equation l|13|l . 
Therefore, surprisingly, the variational potential turns 
out to be self-consistent and hence to provide us with 
the exact baryon of the massive GN model in the large 
N limit. We have been rather explicit here because al- 
though our variational calculation agrees with Ref. [T^ . 
we have come to a different conclusion concerning the 
self-consistency of this ansatz. 

Let us summarize our findings. In the massless limit, 
the shape of the self-consistent scalar potential depends 
only on the filling of the valence level. For complete filling 
{i> — 1), the well degenerates into infinitely far separated 
kink and antikink. As one decreases z/, the baryon shrinks 
and eventually approaches a 1/ cosh^-shape as expected 
from the non- linear Schrodinger equation |l5j| . In the 
massive Gross-Neveu model, we have another parameter 
7 which controls the shape of the scalar potential. If we 
consider only complete filling and turn on the bare mass 
term respectively 7, the baryon also shrinks and even- 
tually becomes non-relativistic. This is physically very 
reasonable. Interestingly, no new shape of S{x) appears 
as we turn on the bare quark mass. It is this feature 
which makes the massive GN model exactly solvable, a 
fact which we had not expected. The only question is 
how the parameter 9 in the general kink-antikink poten- 
tial depends on 1/ and 7. This is answered (implicitly) 
by Eq. (|14|) . Eq. (|15|l then yields the baryon mass as 
a function of fermion number ly and symmetry breaking 
parameter 7. 



III. BARYONS IN THE MASSIVE NJL2 MODEL 
— DERIVATIVE EXPANSION 

In the present section, we consider the NJL2 model 
with continuous chiral symmetry. Unlike in the case of 
the GN model, here we have not been able to guess the 
self-consistent Hartree-Fock potential for arbitrary bare 
quark masses. The only clue we have is the fact that near 
the chiral limit, the sine-Gordon model should somehow 
emerge. In this limit, we know that the relevant length 
scale over which the potential varies is given by l//i, the 
inverse pion mass. This suggests that the theoretical in- 
strument of choice should be the derivative expansion, 
at least in the vicinity of the chiral limit. We will see 
that the sine-Gordon equation is nothing but the lead- 
ing order term of such a derivative expansion. Moreover, 
the calculation of higher order corrections can be done in 
closed analytical form, so that we will get a quantitative 
understanding of the baryons in a reasonable range of 
bare mass parameters, although not for arbitrary mass. 
This is a nice example where one can study in all detail 
how a Skyrme-type description of baryons emerges from 
an underlying fermionic theory, including an unusually 
good control of the corrections. 

As far as the derivative expansion is concerned, we 
refer the reader to the original paper [l^ as well as to 
Ref. ^3 from which we take over the basic formulae 
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and the notation. Nevertheless, we shall try to keep our 
paper self-contained as far as possible. 

To avoid confusion, we should like to draw the at- 
tention to the following difference between this section 
and the preceding one. In the case of the NJL2 model, 
the Skyrme-picture applies only to the case where the 
fermion number is an integer multiple of N, the num- 
ber of flavours (this correponds to = 1 in Sect. II). 
One should not think of filling a positive energy valence 
level, but rather that the winding number of the potential 
induces fermion number by sending energy levels from 
positive to negative energies or vice versa, changing the 
spectral asymmetry. This is discussed in more detail in 
Refs. Only for this case of complete filling of sin- 

gle particle states, the baryon gets massless in the chiral 
limit. One can also contemplate partially filled states in 
the NJL2 model, along the lines of the original work of 
Shei 18], but we shall not consider this possibility here. 
This allows us in the following to define the baryon num- 
ber as an integer, namely fermion number divided by N. 



A. General setup of the derivative expansion 

The central quantity for computing ground state en- 
ergy and baryon number for a system with Hamiltonian 
H is the spectral density 

a{E) ^ TTd{H- E) ^ -ImTri?(£; + ie), (21) 

TT 

where we have introduced the resolvent 



with 



i?(z) = Tr--^^Tr-i^. 



(22) 



The induced fermion number is (—1/2) times the spectral 
asymmetry, 

W = -\j dE<j{E)sgn{E) (23) 

= Im / dE [R{E + ie) + R{-E - ie)] , 

27r Jo 

whereas the ground state energy can be written as 
{H) = I dEEa{E) 

J —00 

1 /"O 

= -Im / AEER{E + ie) (24) 

In a notation similar to Ref. ^3 we decompose H and 
formally into 

H = K + I , H'^ = Hl + V 



(25) 



and expand the resolvent in powers of V ^ 

1 



R{z) = Tr(X + / + z) 



Hi- z^ + V 
Tr(X + I + Z) (Go Y^i-VG^r] (26) 



Go 



(27) 



If one commutes the F's through the Go's by repeatedly 
applying the identity 



GoV^VGo + Ga [V,Hl] Go, 



(28) 



one generates the derivative expansion (the commutator 
[V, Hq] involves derivatives of V). As is well known, this 
method quickly becomes tedious due to the proliferation 
of higher order terms. Therefore we shall use another 
technique based on momentum space below. 



B. Specialization to the massive NJL2 model 

In the large N limit, the Hartree-Fock approach is ad- 
equate. We can therefore choose for H a single particle 
Dirac Hamiltonian with general scalar and pseudoscalar 
local potentials. Without loss of generality, it can be cast 
into the convenient form 

H = -7^19:^ + ^07° 

+ [m + \{x)] e'T'x(:'^)^0g-i7^x(:'=) (29) 

which exhibits the bare mass term and a Hartree-Fock 
potential equivalent to 



with 



(T(a;)7° -|- 7r(x)i7"'" 



a = (m-|-A)cos2x 
TT = — (to + A) sin 2%. 



(30) 



(31) 



We have introduced a mass parameter to which differs 
from the physical fermion mass to. 



m = TO -|- mo 



(32) 



to ensure that A(a;) vanishes asymptotically. The 7- 
matrices will be taken in the representation 7'^ = di , 7^ = 
— i(T2,7^ — f3- Next we identify the operators K and / 
from the preceding subsection as follows. 



H = K + I 



( -id. 


' V 






\ 


id, J 







with 



$ = (to -I- A)e^'^ 4^ Too 



(33) 



(34) 



A natural way of decomposing into Hq and V is 

Hl={-dl + m^) (35) 



V = 



i($')* 



-i$' 



(36) 
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where, using Eq. 134|l . 



I$l 



= 2mA + + 2too('7i + A) cos 2% — 2moTO 



(A' + 2i(m + A)x') e 



(37) 



We shall perform the derivative expansion of the basic 
building blocks 



Tt{K 



z)Go{VGQr 



(38) 



of the resolvent as follows. Consider first the term pro- 
portional to z. In momentum space, 



TvGoiVGo) 



dpdqi dqn-i ^ 



GoipYGoip + qi).. 



2tt 2tt 2tt 

Go{p + qn-l) tr V{qi)V{q2~qi)...Vi~qn-l) (39) 

where now tr is only the Dirac trace. In p-space, the 
potentials vary rapidly as compared to the Green's func- 
tions. We can therefore expand the product of Go's in 
a power series in the g^'s. Once this is done, we trans- 
form the potentials back to coordinate space and carry 
out most of the integrations. The g^'s are replaced by 
derivatives acting on the V' s, 

TrGoiVGoT ^ J '^^ J ^Go{p)^Go{p + qi)... 

Go(p + g„-i)|g^^,(g^+...+a^) 
■ti:V{x,)V{x2)...V{x,,)l^^^. (40) 

In this process, partial integrations are carried out which 
can be justified for the functions A, x which will come 
out at the end (there are no surface terms). We have 
suppressed the Taylor expansion in Eq. H4U|) in order to 
keep the structure of the formula transparent. The other 
two terms in Eq. I|38|l can be handled similarly with the 
result 

Tr/Go(FGo)" = / / ^Go(p)Go(p + ft)... 

Go(p + <?„)|,,^,(a,+...+9,) 

•tr I{x,)V{x2)V(xs)...V{xn+i)l^^^ 

TrKGoiVGor = J "^"^ J ^pGo{p?Go{p + qi)... 

Ga{p + qn^l)\q^=,^a,+,„+d^) 

■tr <j3V{xi)V{x2)...V{xn)l^^, . m 

If we apply the above formalism to Hartree-Fock, we must 
of course add the double counting correction to the en- 
ergy. 



E^(H) 



dx 



(m + xy 



(42) 



C. Energy 

Only the odd part of the resolvent R(z) in Eq. H26|l con- 
tributes. In the expansion in powers of V, the terms of 



order and play a special role. They contain diver- 
gencies which require regularization and renormalization, 
whereas all higher order terms are finite. Besides, they 
do not necessitate the derivative expansion (there are no 
qi to expand in). We therefore show this calculation in 
some detail first and then discuss the higher order terms. 
Terms of order V'^ 

Here we are considering the vacuum energy. Although 
this will be subtracted from the baryon energy, it is useful 
to calculate it in order to rederive the renormalization 
condition needed later on. The resolvent is given by 



R{z) = dx 

J J-A/2 



dp 



2z 



2TT (p2-HTO2) 



(43) 



with an UV cut-off A/2. Eqs. ^ and (gH then yield 
the energy density 



dp r-^ . 

— VP + + 



A/2 



27r 



2Ng^ 



= -8^-4^+2^1^^A+S]V^- ^44) 

Minimizing with respect to the dynamical fermion mass 
TO, we recover the well-known gap equation 



(to — Too) to ^ A 

Ng'^ TT TO 



(45) 



Upon using this condition, the vacuum energy density 
becomes (up to an irrelevant quadratic divergence) 



TO 

47r 



totoq a 

In — , 

27r TO 



(46) 



As pointed out above, toq /Ng"^ or equivalently Too In A/to 
is a physical quantity. Here, it is more natural to express 
it in terms of the pion mass /i (see ^), 



TTTOo 



1 



: arctan ■ 



1 



V 



For later convenience, we introduce the function 
4 y 



(47) 



Fiy) 



y 



arctan ■ 



6-^ 30'^ 140^ 



(48) 



which allows us to express the renormalized vacuum en- 
ergy density in the following way, 



^vac = - ^F{^,/m). (49) 



2 ^2 



This already shows that an expansion for small quark 
masses should be thought of as an expansion in the ratio 
of pion mass to (dynamical) quark mass in the NJL2 
model. The bare quark mass TOq goes to in the limit 
A 00 and cannot appear in any physical quantity. 
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Terms of order 

Now we have to compute 



dx 



dx ■ 



dEElmTrzGoVGo 

o 

(m + A)^ — 
1 /-^/^ dp " 



■ImzGlip) trV{x) 



dEE 
{fa + A)2 — 



(50) 



E. Results 

On the basis of the effective sine-Gordon theory for 
the NJL2 model, the mass of the baryon with fermion 
number N is expected to approach 2/i/7r for /i — > 0. Our 
goal is to compute corrections to this result up to or- 
der /i®/m®. This requires a substantial number of terms 
in the derivative expansion and correspondingly involved 
algebraic expressions. To simplify the notation, let us set 
m = 1 from now on. Then, the energy density in the rel- 
evant order of the derivative expansion computed along 
the lines of Sect. Ill C has the following form, 



The integrations over E and p can easily be carried out. 
Owing to the gap equation, all divergencies cancel out 
and the unphysical quantities {Ng^ , mo. A) can again be 
eliminated in favour of /i and m with the result 



dx 



(51) 



-\ (l-cos2x) + — 
m / 2m? 



Once again we have dropped all terms in which the bare 
mass mo is not multiplied by a factor ln(A/m). 
Terms of order , n > 1 

All higher order terms in the expansion (|26|l are free of 
UV divergencies. Therefore, we can set the bare quark 
mass mo — 0, thereby greatly simplifying our potential 
V , and calculate the traces and integrals mechanically 
with computer algebra (we used Maple). The result for 
the energy density coming from these higher order terms 
is a polynomial in A, x and their derivatives. Detailed 
results will be given in Sect. Ill E. 



2-k£ 



{l + X) (cos2x-1)--A2 



+(x')'-^(x")' + ^(x"')' 
^:^(x")^ + A^ + ^(A')^ 



140 



(X 



IV\2 



-A* 



120 



(A' 



//\2 



ll\2 



^a(x"')^ + JaxV^ 



A^(X")' (53) 



We have subtracted the vacuum contribution and sim- 
plified the result as much as possible with the help of 
partial integrations (only the energy Jda;£ is uniquely 
defined in this approach). To appreciate the complex- 
ity of expression (|53|l . one should compare it with the 
leading order terms only. 



2n£ = (x') 



l\2 



-(cos2x-l). 



(54) 



D. Baryon number 

The leading order contribution comes from the 
TtIGqVGq term in Eq. (|26() and can easily be shown 
to have the expected topological form 



{N) 



dx^^-[x{^)-x{~^)], (52) 

TT TT 



with the well-known identification of winding number of 
the chiral phase with baryon number. Since we were not 
sure to which extent the general topological arguments 
carry over to a system with explicit symmetry breaking, 
we have computed the next three orders in the derivative 
expansion, using Maple. We found that the result 1)52(1 
does not get any correction whatsoever and therefore pre- 
sumably holds to all orders. Although many terms are 
produced in the integrand (the baryon density) by our 
algorithm, they can be nicely combined into total deriva- 
tives of functions which vanish at infinity and therefore 
do not affect the lowest order topological result (|52|l . 



which reproduce exactly the sine-Gordon theory and rep- 
resent the state of the art prior to this work. We then 
vary the energy with respect to x and A. In order to solve 
the resulting differential equations, we expand -F(/i) and 
assume the following Taylor series for x and A, 



X ~ Xo + M^Xi + M^X2 + M^X3 
A « /i'Ai+//A2+/i''A3. 



(55) 



The coefficients Xn and A„ in turn are taken to depend 
only on ^ = ^x, so that each derivative with respect 
to X increases the power of /i by one. All of these as- 
sumptions can be justified a posteriori by showing that 
they are the simplest ones which lead to a consistent 
approximate solution of the differential equations. This 
procedure yields the following set of inhomogeneous dif- 
ferential equations for Xn and algebraic equations for An 
(where now ' = d/d^). 
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Ai 

Xi - X\ cos 2x0 
A2 



X3 - X3 COS 2x0 



; sin 2x0 

- (cos 2x0 - 1) 

; (Ai + i)sin2xo-^xr 



= ^ (cos 2x0-1) 



j(Xo) 



X2-X2cos2xo = 



2X1 sin 2x0 

4^^ 2^1 + 12^1 
X? + + ^Ai + 1) sin 2x0 

-1X1 + ^Xi^ cos 2x0 - ^xS'^ 

'3A1X0 + gAlXo + 3A1X0 - gXl 



+ A 



A, = 



-AlA; 

1 



// // 

■XoXi 



20-Yo^oXo -3o(Xo) - 4^2 + Y2^2 + gAi 

-(A'i)2 - IaiA'/ + — A^ - -Ai + -Ai(xf,')' 
\) 6 ^ 120 ^ 24 2 

^ ^ 2.1 1. \ „ 

I gX2 + A1X2 + A2XI - gXl + ggXl + g 1 X 1 I COS2X0 



1 

140 
2 



V 

1 ,in^O 



280 
yiu 4 



--Xo(Xo) 



-j^glXoJ Xo +2-^1X0 +3^1X0 

— AiXo + 5A1X0 + Y^Ai Xo - gX2 + 3-^2X0 + 3-^1X1 

^ 3 -^1 3^^2X0 + 3 AiXi + 3^2X0 + 3^^2X0 

^-^1X1 - ^Xi 



-^Alx 



^Aia;'x^'-2Aia;x[,"-^ 



A?xr - (A'i)^X^' 



(56) 
(57) 
(58) 

(59) 



(60) 



(61) 



(62) 



In spite of their frightening appearance, it is not too hard 
to solve these equations recursively in the order in which 
they have been written down here. Our analytical results 
for baryon number 1 are (sech^ — ^/ cosh^) 



1 



X3 



Xo 
Ai 

Xi 
A2 

X2 
A3 



— 2 arctan e^ 

= — sech^f 
2 ^ 

1 2 
= - sinh ^ sech ^ 
8 

= — - (sech^^ — sech^f) 

= - -^r- sinh C (23 sech^C - 122 sech''^) 
1152 ^ ^ 

= -TTTTT (252 sech^C - 1225 sech"*^ + 1061 sech^C) 
1440 ^ ^ 



691200 

-90456 sech^C) 



sinh,^ (2621 sech^^ 

6, 



108092 sech^C 

(63) 



In the course of this calculation, we had to decide what 
to do with the homogeneous solution of the differential 
equations for Xn- Physically, they reflect the fact that 
there is a flat direction in function space due to the break- 
ing of translational invariance by the baryon. We have 
made the choice that x is odd under ^ — > — ^, a require- 
ment which fixes the position of the baryon in space and 
leads to a unique solution of the differential equations. 

We see that the leading order results (keeping only xo) 
agree with the sine-Gordon theory, whereas the higher or- 
der corrections yield systematic corrections to it. Insert- 
ing Eqs. H63|) into Eqs. H55|l . the scalar and pseudoscalar 
potentials can be written as a power series in /i^, i.e.. 
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the ratio of pion mass to physical quark mass (remember 
that we have chosen units such that m — 1), with smooth 
coefficient functions depending only on fix. 

Inserting the results into the expression for the 
energy density H53|l and integrating over dx, we finally 
get the following chiral expansion for the baryon mass in 
the NJL2 model, 

(64) 

Due to the values of the coefficients, the corrections are 
numerically small even for /i/m ~ 1. We see no difficulty 
in principle to push this calculation to higher orders, but 
the number of terms appearing in intermediate steps of 
the calculation increases rapidly. It then becomes in- 
creasingly difficult for Maple to handle and simplify the 
lengthy expressions. 

IV. SUMMARY AND CONCLUSIONS 

The purpose of this work was to begin a systematic 
study of massive Gross-Neveu models in 1+1 dimensions. 
As a warm-up for a comprehensive study of the phase 
diagram of these models, we have reconsidered the issue 
of baryons. Two different strategies turned out to be 
successful for the massive discrete chiral GN model and 
the massive continuous chiral NJL2 model, respectively. 

In the case of the GN model, we simply guessed the 
scalar potential for a Hartree-Fock calculation and proved 
its self-consistency. Actually, the shape of the potential 
is the same as in the massless limit, but for a differ- 
ent fermion number. This was crucial for being able to 
carry through the calculation, since potentials for which 
one can solve the Dirac equation analytically are ex- 
tremely rare. The same calculation in a somewhat differ- 
ent framework has already been done some time ago 12]; 
however, the decisive fact that the result is self-consistent 
was apparently missed by the authors. This baryon solu- 
tion should be a good starting point for studying hot and 
dense matter in the GN model, generalizing our recent 
work in the chiral limit Q to finite bare quark masses. 

Baryons in the NJL2 model are a totally different story: 
In the chiral limit, they become massless, baryon number 
has a topological meaning, and a Skyrme-type physical 
picture is adequate. Here we applied a derivative expan- 
sion technique which allows us to bypass the cumbersome 
Hartree-Fock procedure. Without explicitly solving the 
Dirac equation, we have nevertheless obtained quantita- 



tive results encoded in the first three correction terms 
beyond the sine-Gordon limit. In effect, this procedure 
can be regarded as a kind of bosonization where a chiral 
angle field x £^nd a radial field A carry all the dynami- 
cal information. The emergence of the Skyrme picture is 
put on very solid grounds in this case. The small parame- 
ter which governs the expansion of potentials and baryon 
masses is the ratio of pion mass to dynamical quark mass, 
whereas the bare coupling constant and the bare quark 
mass disappear in the process of renormalization. 

Unfortunately, it is unlikely that this method carries 
over to higher dimensions. The fact that the size of 
the baryon increases like the inverse pion mass in the 
NJL2 model was instrumental here, but is not expected 
to hold in more than one space dimension. Nevertheless, 
we hope that our results are of some use to test algo- 
rithms or study questions related to chiral perturbation 
theory with baryons. 

Finally, let us comment on the relationship between the 
GN model and quasi-one-dimensional condensed matter 
systems. As is well known and has recently been re- 
emphasized by us [^] , the GN model can serve to describe 
systems like the Peierls-Frohlich model, one-dimensional 
superconductors, or conducting polymers. Does the step 
from the massless to the massive GN model have any 
analogy as well? The answer is yes — it is the transition 
from systems with a two-fold degenerate ground state to 
non-degenerate systems. The prime example is perhaps 
the step from trans-polyacetylene to cis-polyacetylene . 
Following the work of Brazovskii and Kirova |l9j , a num- 
ber of studies have been devoted to the existence of 
polarons, bipolarons and excitons in such doped, non- 
degenerate polymers (doping changes the number of elec- 
trons in the valence band), see the reprint volume [20j| . 
We can compare these studies to our results in Sect. II for 
N = 2 (there is no flavour in the condensed matter case, 
only two spin states). Thus, the bipolaron corresponds 
to our V = 1 case, the polaron to v = 1/2, and indeed 
for these values the mathematics of the condensed mat- 
ter continuum models is identical to our calculation. In 
the polymer case, one also considers states with one elec- 
tron in the negative and positive energy valence levels, an 
exciton. As a matter of fact, our calculation can be re- 
peated for arbitrary occupation fractions i>± for these two 
discrete levels, preserving self-consistency. Such states 
would correspond to heavier, excited baryons in the GN 
model, with fermion number N{vj^ + — 1) and an ad- 
mixture of anti-quarks. If the baryon number vanishes 
(I'-i- — 1 — !/_), one could perhaps talk about "baryo- 
nium" as the analogue of an exciton. 
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